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NOMENCLATURE 


Symbol 

d 

G 

LWC 

M w 

T 

a 

T f 

T 

o 

V 

V 

e 

W f 

X 

Greek 

Symbol 

o 

S 


Definition 

Water droplet diameter 

Glycol mass fraction in glycol-water solution 

Liquid water content in atmosphere, gm/m 

Mass rate of water droplet impact per unit area on wing 
surface 

Ambient freestream temperature 

Freezing temperature of glycol-water solution 

Stagnation temperature 

True airspeed 

Equivalent airspeed 

Rate of flow of glycol solution per unit area 

Glycol mass fraction in glycol-water solution pumped 
through porous leading edge 


Angle of attack 
Local catch efficiency 


Subscripts 

L Lower panel 
U Upper panel 
max Maximum 
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INTRODUCTION 


The concept of pumping a glycol-water solution through a porous 
leading edge skin to achieve ice protection is not new. Operational 
systems using this concept have been employed on European airplanes 
for many years. However, no U.S. manufacturers have used liquid anti- 
ice protection for wing and tail leading edges up to this time. Re- 
cently there has been a surge of interest in the U.S. in this concept 
because of the advantages it offers in comparison with pneumatic boots 
and hot air bleed systems. 

Several icing tunnel tests with liquid ice protection systems have 
been conducted in the NASA Lewis Icing Research Tunnel during the past 
few years to add to a relatively meager data base. The purpose of this 
report is to present a proposed method of analytically predicting the 
minimum fluid flow rate required to provide anti-ice protection with a 
porous leading edge system on a wing under a given set of flight condi- 
tions. Results of the proposed method are compared with the actual re- 
sults of an icing test of a real wing section in the NASA Lewis Icing 
Research Tunnel. This work was conducted under NASA Grant NAG 3-71. 

SYSTEM DESCRIPTION 

A fluid ice protection system consists of a porous leading edge 
skin panel attached to the leading edge of a wing and a pump that dis- 
tributes a glycol-water solution from a reservoir to the leading edge 
panel through nylon tubing. The fluid flows through the panel onto the 
surface of the wing, providing either an anti-icing capability by dis- 
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solving supercooled water droplets and preventing the formation of ice, 
or a de-icing capability by chemically breaking the bond of established 
ice. A significant feature of the system is that protection is obtained 
aft of the panel by the flow of the fluid along the chord to the trail- 
ing edge, thus preventing the formation of ice anywhere aft of the ac- 
tive leading edge. 

The porous panel material most commonly used consists of two or 
three layers of stainless steel wire cloth that are rolled, sintered, 
and finish rolled to proper thickness. Recent development programs 
have also produced porous panels made of laser-drilled titanluir sheet, 
and various composite materials. 

A typical cross section of a porous leading edge panel installa- 
tion is shown in Figure 1. The edges of the active portion of the 
panel must be placed such that extreme positions of the stagnation 
point for which icing protection is required are not too close to the 
edge to prevent fluid from being distributed on both the upper and 
lower surfaces of the wing. 


PREDICTION METHOD 


Anti-ice protection is obtained by providing a glvcol-water solution 

on the leading edge of the wing that mixes with the atmospheric water 
droplets as they impact on the leading edge. To provide anti-ice pro- 
tection, the resulting solution must have a glycol mass fraction that 
is high enough to prevent any freezing at the leading edge or on the 
wing upper and lower surfaces as the fluid flows aft along the wing. 
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Figure 1: Cross Section of a Porous Panel Installed 

on a Wing Leading Edge. 
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original page is 

OF POOR QUALITY 

The water droplets do not Impact as a uniform mass rate on the lead- 
ing edge. The mass rate per unit area tends to be highest at or near the 
stagnation point and decreases in a chordwise direction on either side of 
the stagnation point. The distribution of water mass rate on the leading 
edge is described by the nondimensional catch efficiency parameter 6. 

The actual local water mass rate of impact per unit area is given as 

M w = B(LWC) (V) . 

A typical distribution of B is shown in Figure 2. The trajectories of 
the water droplets relative to the flow streamlines determine the B 
distribution. Thus S is a function of the airfoil shape and size, the 
remote airspeed, the air density, and water droplet diameter. Given S, 
the local and total mass rate will be directly proportional to liquid 
water content . 

Until recently, the only methods of predicting S were empirical 
and semi-empirical. However, the development of improved methods of 
computational aerodynamics has resulted in several different computer 
codes that predict 8 distributions with relatively good accuracy (Ref- 
erences 1 and 2). 

In contrast to the 6 distribution, the glycol-water fluid pumped 
through the porous leading edge is distributed at a uniform rate over 
the surface of the porous skin as shown in Figure 3. 

The basic assumption of the proposed prediction method is that the 
minimum glycol flow rate that will still achieve anti-icing (no ice 
accumulated on the leading edge), results in a freezing temperature of 
the glycol-water mixture, at the point of maximum water catch rate, 
that is equal to the local air temperature. At this point on the wing 
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(near or at the stagnation point) , the glycol concentration is mini- 
mized; and at the anti-ice threshold the mixture will just begin to 
exhibit the formation of small pieces of ice. As the glycol flow rate 
is decreased, the extent of ice formation gradually increases until 
continuous bars are formed spanwise on the leading edge before being 
swept away every few minutes. This latter mode is called natural 
deicing. 

Obviously, to determine the appropriate glycol rate, the local air 
temperature must be known. This temperature will be between stagnation 
temperature and ambient atmospheric temperature. To be conservative, 
one would choose ambient temperature because it will require a lower 
freezing temperature and higher flow rate. However, if a constant flow 
rate is chosen based on the most severe icing conditions anticipated, 
it is probably sufficient to use stagnation temperature as the local 
glycol-water freezing temperature. In this case, the system will per- 
form somewhere between an anti-ice mode and natural deice mode during 
the most severe conditions. Reference 3 shows that the flow rate re- 
quired for natural deicing is only 25% to 50% of that required for 
anti-icing. At less severe icing conditions, the system will have a 
flow rate equal to or greater than that i quired for anti-icing. As 
an appropriate compromise, it is suggested that the mean temperature 
between ambient and stagnation be used. 

The method is utilized by following these logical steps. 

1. Find 8 as a function of 
max 

- airfoil shape 

- airspeed 
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- air density 

- droplet diaEkier. 

Any reasonably accurate computer code, such as In Refers aces 

1 and 2, can be used for this step. 

2. Calculate the water catch rate N by the formula 

w 7 


N 


w 


6 


(LWC) <J. 


3. Determine the glycol mass fraction G required to produce a 
solution with a freezing temperature equal to the average be- 
tween ambient static temperature and stagnation temperature, 
using the graph shown in Figure 4. 

4. Calculate the fluid flow required to achieve the glycol mass 
fraction G, given a water catch rate M w , bv the equation 

GM 

w 

"f X - G 


where X is the initial glycol mass fraction. For most fluid 
systems X is approximately 0.8. 


PREDICTED AND ACTUAL FLOW RATES 

To test the validity of this method, it was applied to an airfoil 
for which anti-ice fluid rates were determined by tests conducted in 
the NASA Lewis Icing Research Tunnel and reported in Reference 3. 

Those tests utilized an actual airplane wing section. The orig- 
inal wing tapered from a NA r A 642A215 airfoil at the root (WS 0) to a 
NACA 64jA412 at the tip (WS 216). The wing incorporated increased 
thickness on the forward 30% of the upper surface, a modification pro- 
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Figure 4: Freezing Temperature of a Monoethylene Glycol-Water Solution. 
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posed by R. Hicks of NASA Ames Research Center (Reference 4). The 
centerline of the tunnel was at VS 58 of the original wing, where the 
chord is 63.25 in. The predictive method was applied at this station. 

The porous panel mounted on the leading edge consisted of three 
independently controlled sections, with WS 58 in the center of the 
center section. To obtain the anti-ice flow rates, the upper and 
lower sections were used simultaneously during each run to establish 
independent flow rate values while the center section was used to ob- 
tain minimus flow rates for natural deicing. 

The method of determining the glycol flow rate corresponding with 
the anti- ice threshold was as follows. At a given flight condition, 
the flow rate was set to be well above the anti-ice threshold. The 
flow rate was then reduced in steps, allowing about 30 seconds for the 
system to stabilize at each point, until small flecks of ice began to 
appear on the leading edge in the vicinity of the stagnation point. 

At the anti- ice threshold, the small ice flecks, ranging up to about 
3 mm in diameter, would form and then be swept downstream in only a 
few seconds. A glycol flow rate lower than the threshold value would 
allow the ice flecks to persist and to grow gradually into larger patches 
before being shed from the wing. Complete test details may be found 
in Reference 3. 

Six test cases were chosen for analysis. Test conditions and flow 
rate data are presented in Table 1. Also shown are the computed values 
of 8 mny . The computer code developed at Ohio State University (Ref- 
erence 2) was utilized. 
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Table 1: Test Data and Predicted Maxinun Catch Rate 

for an Anti-ice Porous Leading Edge System 


Case 

V e 

kt 

LUC 

gm/m 3 

d 

urn 

T o 

°F 

a 

deg 

W f 

U 

ml/cmmin 

u 

f L 

ml/cm 2 mln 

B max 

1 

96 

1.50 

15 

25 

7.8 

.0255 

.0193 

.208 

II 

96 

2.40 

20 

25 

7.8 

.0436 

.0340 

.314 

III 

96 

1.16 

15 

5 

1.2 

.0420 

.0330 

.219 

IV 

175 

.65 

11 

25 

1.2 

.0240 

.0250 

.259 

V 

175 

.80 

15 

25 

1.2 

.0380 

.0350 

.329 

VI 

175 

1.16 

15 

25 

1.2 

.0540 

.0430 

.329 


To check the validity of the proposed method, the freezing temper- 
ature of the glycol-water mixture at the stagnation point, given S and 
the glycol flow rate, was determined by the following method. 

The water catch rate M was calculated as 

w 

M w - B(LWC)V ; 


then 


G 


XW, 


W* + M 
f w 


For these tests, X * 0.8. Knowing G, one can determine the freezing 
temperature of the glycol-water mixture using the data shown in Figure 4 
To determine the sensitivity of the freezing temperature to 6, cal- 
culations were made for a range of B above and below the computed value 
of S njax for each case. The sensitivity of freezing temperature to the 
glycol flow rate was determined by varying above and below the experi 
mentally observed values of at the anti-icing threshold by 205. 
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Results of the calculations are presented in Figure 5 for the six 
test cases. 

Curves are presented for the upper and lover panels. The center 

panel, where 5 was calculated, would be expected to produce a curve 

between the upper and lower panel. The consistently lower predicted 

freezing temperature of the upper panel fluid at the stagnation point 

is caused by the fact that the calculated 8 underpredicts 8 for 

aax sax 

the sharper leading edge of the upper panel of the tapered wing. The 
reverse is true for the lower panel. 

In every case except for Case I, the midpoint between the upper 
and lover panel curves predicts a freezing temperature at the point 
of £ that lies between the ambient and total air temperature , a 
result to be expected. This calculated freezing temperature is based 
on the predicted value of S, the observed glycol flow rate, the wind 
tunnel test conditions, and the properties of a glycol-water solution. 

The sensitivity of the predicted freezing temperature to errors 
in the value of S and to variations in the glycol flow rate are illus- 
trated for each case. 

It is now useful to compare actual values of W ^ with those pre- 
dicted by the method presented herein. Results are shown in Table 2. 
Note that the glycol flow rate is converted from volume flow rate to 
mass flow rate by the specific gravity of the original fluid, which 
is 1.1. 
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v e * 96 

CASE // fc- 25 

LWC * 2.4 
*C- 7.8 




Figure 5: Effect of 6 on Local Fluid Freezing Temperature, and 

Effect of Flow Rate on Local Freezing Temperature at 

Predicted 6 . Case II. 

max 
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STfOOR ***•*" 


v 9 = 175 

CASE IV r o = 25 

LWC - .65 
<< - LEO 




Figure 5: Effect of 8 on Local Fluid Freezing Temperature, and 

Effect of Flow Rate on Local Freezing Temperature at 

Predicted 6 . Case IV. 

max 
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Table 2: Comparison of Predicted with Actual Test Flow Sates 

at Anti-icing Threshold 


Case 

8 

max 

M w 

gm/cm^min 

T + T 
a o 

W f 

predicted 
gm/ cm'-mic 

W + w 
U L 

Prediction 

Accuracy 

% 

2 

deg F 

2 

go/cm 2 min 

1 

.208 

.0837 

23.9 

.0179 

.0246 

-27.2 

11 

.314 

.202 

23.9 

.0432 

.0427 

+ 1.2 

III 

.219 

.0666 

3.9 

.0415 

.0413 

+ .5 

I\ 

.259 

.0878 

21.5 

.0228 

.0269 

-15.2 

V 

.329 

.137 

21.5 

.0356 

.0401 

-11.2 

VI 

.329 

.199 

21.5 

.0517 

.0533 

- 3.0 


These results show that the method of prediction of anti-ice thresh- 
old flow rates presented herein predicts flow rates within an average 
error of less than 10 percent of the experimentally determined flow 
rates. This is believed to be excellent, considering the accuracy of 
the S prediction and the fact that the anti-icing threshold character- 
istics tended to persist over a relatively wide range of values of glycol 
flow rate, making it difficult to obtain fine resolution of the anti- 
icing threshold. 

Therefore, this method appears to provide a reasonably accurate 
determination of the flow rate required to assure continuous anti-icing 
performance at a given icing flight condition. For less severe icing, 
an excess of anti- ice fluid is available at the leading edge. For 
more severe icing, there will be a gradual transition to a natural 
deicing mode which still provides adequate protection against ice 
buildup in flight. 
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It is recommended that additional testing be done on different 
airfoils under a variety of test conditions to verify the generality 
of the method presented in this report. 
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